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Multi-terminal Aharonov-Bohm (AB) rings are ideal building blocks for quantum networks (QNs) 
thanks to their ability to map input states into controlled coherent superpositions of output states. 
We report on experiments performed on three-terminal GaAs/Al^Gai-xAs AB devices and compare 
our results with a scattering-matrix model including Lorentz forces and decoherence. Our devices 
were studied as a function of external magnetic field (B) and gate voltage at temperatures down 
to 350 mK. The total output current from two terminals while applying a small bias to the third 
lead was found to be symmetric with respect to B with AB oscillations showing abrupt phase 
jumps between and ty at different values of gate voltage and at low magnetic fields, reminiscent 
of the phase-rigidity constraint due to Onsager-Casimir relations. Individual outputs show quasi- 
linear dependence of the oscillation phase on the external electric field. We emphasize that a 
simple scattering-matrix approach can not model the observed behavior and propose an improved 
description that can fully describe the observed phenomena. Furthermore, we shall show that 
our model can be successfully exploited to determine the range of experimental parameters that 
guarantee a minimum oscillation visibility, given the geometry and coherence length of a QN. 

PACS numbers: 05.60. Gg, 73.23.-b, 73.63.-b, 85.35.Ds 



I. INTRODUCTION 

A quantum network (QN) can be implemented by 
means of a set of nodes connected by electron waveguides 
that allow the coherent control of electron wavefunctions 
by means of external magnetic or electric fields. One of 
the simplest QNs based on coherent electrons is a three- 
terminal Aharonov-Bohm (AB) ring with one input and 
two output channels^. When the linear superposition 
of the two output states can be tuned by means of an 
external field, this device realizes a qubit. The imple- 
mentation of single-qubit logic functions was proposed 
based on appropriately tailored multiterminal-terminal 
rings 1 2 , making these QNs promising building blocks for 
quantum-computation architectures. Devices were the- 
oretically analyzed by means of scattering-matrix ap- 
proaches in the Landauer-Biittiker frameworlP^, free- 
electron-like node equation d 1 * 2 * 5 1 , and displaced Gaussian 
wavefunctions 6 . The effect of Lorentz forces in three- 
terminal rings was theoretically studied in Ref. 6 by solv- 
ing the time-dependent Schrodinger equation. Despite 
this intense research activity, to the best of our knowledge 
these approaches have not yet been compared against ex- 
perimental results on real devices. 

In this work we shall extend existing theories based 
on scattering-matrix approaches by including the effect 
of decoherence and classical (Lorentz) forces in the de- 
scription of the system and carry out this comparison 
by studying the low-temperature the coherent-transport 
properties of three-terminal AB rings. Our analysis 
shows that the inclusion of these effects is necessary to 
fully understand the details of the observed behavior. 
Additionally, we shall identify a set of criteria that a 
device must satisfy to implement a functional QN and 
provide the operational range that guarantees a chosen 



value for oscillation visibility for a given set of device 
parameters and electronic coherence length. 

Our devices were fabricated starting from a 
GaAs/Al^Gai-^As heterostructure containing a 
high- mobility two-dimensional electron gas (2DEG). We 
observed AB oscillations in both output channels and 
in the total output of the system as a function of an 
external perpendicular magnetic field B. Oscillations 
show a non-trivial phenomenology as a function of gate 
voltage V g \ they evol ve in a continuous way from a 
phase-rigid regime^^, where the phase of the oscilla- 
tions as a function of B shows abrupt changes from to 
7r when an external electric field E is applied, to a linear 
dependence of the oscillation phase with E. 



II. EXPERIMENTAL RESULTS 

Three-terminal AB rings were fabricated from a two- 
dimensional electron gas (2DEG) confined 90 nm below 
the surface of a modulation-doped GaAs/Al^Gai-^As 
heterostructure. At a temperature T — 4 K the un- 
patterned 2DEG density and mobility were found to be 
2.1 x 10 11 cm -2 and 1.7 x 10 6 cm 2 /Vs, respectively. 

The ring geometry was defined by shallow plasma etch- 
ing. The same processing step realized a set of lateral 
gates (labeled Gi through G5) that provide control over 
the electron density in each ring arm. A SEM image in 
artificial colors of one of our devices is shown in Fig. lb. 
The ring external (internal) radius is 220 nm (90 nm). 
Standard Ni/AuGe/Ni/Au (5 nm/180 nm/5 nm/100 nm) 
n-type Ohmic contacts (not shown in the figure) were 
fabricated to allow electrical access to the 2DEG. 

The AB interference pattern can be tuned by exploit- 
ing the magnetic and/or electric AB effect 1 ^ by means of 
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a perpendicular magnetic field or by changing the volt- 
ages Vgi, Vgi, and Vg3 applied to gates Gi, G2, and G3 
respectively. 

Measurements were performed in a 3 He cryostat at 
350 mK in a three-terminal configuration, where the in- 
jection contact was biased with an ac excitation signal 
V ex = 30 /iV at frequency of 170 Hz and currents I2 
and Is were measured respectively at leads L2 and L3 by 
means of current preamplifiers and phase-sensitive lock- 
in techniques. Blocking capacitors were present at the 
inputs of the preamplifiers to remove any unwanted dc 
component of the bias. In the following we shall refer to 
the ratio gi = Ii/V ex as the conductance of the output 
lead Li, with i = 2,3. Three nominally identical devices 
were investigated in detail and found to display the same 
behavior. In the following we shall focus on one of them, 
discussing when relevant the differences or the similari- 
ties observed in the others. 
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FIG. 1: (Color) (a) AB oscillations measured at a fixed gate 
voltages (Vgi = . . . = Vgb = 0.12V) at 350 mK. The red 
and green lines represent the two output conductance, #2 and 
#3 , while the black line shows the total conductance #2 + #3 • 
(b) SEM image of the ring in artificial colors. The gray area 
corresponds to the 130-nm deep etching. The yellow areas 
represent the five gates (Gi through G5) while the central 
red region is the three-terminal AB ring with one input (Li) 
and two output (L2 and L3) leads. The segment in the figure 
represents a length of 1 /mi. (c) Fourier transform of the 
derivative with respect to B of the total output current shown 
in Panel a. 

Figure la shows #2, 93, and the total output conduc- 
tance g t = #2 + 93 as a function of the perpendicular 
magnetic field at fixed gate voltages. AB oscillations are 
visible in both g^ and gs and in their sum with a maxi- 
mum visibility of approximately 0.2 indicating coherent 
transport across the ring. The period of the oscillations, 
calculated by Fourier transforming the data (Fig. lc), 
is ~ 52.6 mT, corresponding to h/e AB oscillations^, 
where e is the electron charge, and h is the Plank con- 
stant, for a ring with an effective radius of 160 nm, in 
good agreement with the sample geometry. Similar val- 



ues (155 nm and 170 nm) were found for the other devices 
studied. The oscillating part of #2 (#3) is superimposed 
on a background that increases (decreases) as a function 
of B and its amplitude decreases at higher (in modulus) 
magnetic fields. This behavior originates from an im- 
balance in the branching probability due to the Lorentz 
force^ and leads to a suppression of the visibility of the 
oscillations that may limit the operation of a QN and 
should be taken into account in the design of the device. 

The total output conductance was found to be sym- 
metric with respect to the magnetic field as shown by 
the black curve in Fig. la. The second derivative of g t 
with respect to B is reported in Fig. 2 as a function of B 
and Vg2- It does display the same symmetry in the entire 
range of gate voltages explored. In the region between 
-0.1 T and 0.1 T of Fig. 2 abrupt jumps of the oscilla- 
tion phase from to 7r can be seen at Vgi ~ 0.17 V and 
Vgi ~ 0.125 V, reminiscent of the phase-rigidity phe- 
nomena observed in two-lead, closed ring d 9 * 11 *. As the 
modulus of the magnetic field approaches 0.2 T, these 
phase jumps become smoother, and evolve towards an 
almost continuous shift of the phase with gate bias. It is 
interesting to note that, within this continuous-evolution 
range, the oscillations at B > present a phase shift de- 
pendence on the gate voltage that has an opposite sign 
with respect to those at B < 0: the oscillation maxima, 
shown as red areas in Fig. 2, drift to higher (in modu- 
lus) magnetic fields as the gate voltage is increased. We 
stress that this is not due to a change in the frequency of 
the oscillations, in fact the distance between neighboring 
maxima remains approximately constant. 
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FIG. 2: Second-order derivative of the total conductance g<2 + 
gs with respect to the magnetic field plotted as a function of 
B and Vgi- 

A remarkably different behavior was observed in the 
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evolution of the individual outputs as a function of B 
and Vq2, shown in Fig. 3. In this case, the phase of 
the oscillations of g^ evolves almost linearly with Vq2 in 
the entire range of magnetic fields and gate voltages ex- 
plored, gs shows a similar behavior, but with an opposite 
dependence of the phase evolution on Vgi> The same be- 
havior was observed in all the devices studied: the total 
output conductance displayed a transition from a phase- 
rigid regime to a linear regime at \B\ & 0.2 T, while the 
individual outputs displayed a linear evolution pattern in 
all the range of magnetic fields explored. 
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FIG. 3: Contour plots of the second-order derivatives of #2 +#3 
(panel a), g^ (panel b), and gz (panel c) with respect to the 
magnetic field as a function of B and Vg2- 



III. THEORY AND SIMULATIONS 

Here we shall propose a novel scattering-matrix ap- 
proach that is able to explain in detail the observations 
reported in the previous section and will compare our 
numerical results with the experimental data. Following 
Refs. 3,4, electronic transport can be described in terms 
of transmission and reflection coefficients of a N x M 
scattering matrix S that in the most general case links 
M input electron amplitudes (u) with N output ones (u'). 
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In the case of our three-lead devices, the propagation 
of an electron with a given energy is modeled by a 3 x 3 
complex matrix, S, whose coefficients can be calculated 
given the scattering matrices of the three identical blocks 
shown in Fig. 4a that compose the ring. Each individual 
block (Fig. 4b) is composed of a scattering center (black 



triangle in the figure) coupled to three leads. One of the 
leads represents one of the three ring arms (black square) 
where electron phase evolution is driven by wavefunc- 
tion propagation and electric and magnetic AB effects. 
Lorentz forces are modeled by appropriately changing 
the branching probability of the scattering center as a 
function of magnetic field, while decoherence effects are 
introduced in our scheme by decreasing the amplitude of 
the coherent wavefunction of the transmitted electron. 




FIG. 4: (a) Complete scheme of the three-terminal AB ring, 
(b) Representation of a building block describing the coupling 
between one input or output lead and the ring. 

The scattering center is represented by a 3 x 3 scatter- 
ing matrix. Current conservation requires this matrix to 
be unitary. We also restrict ourselves to the case of real 
matrix elements, which is equivalent to assume that no 
phase shifts are introduced by the scattering center. Un- 
der this assumptions the matrix is orthogonal (0(3)) and 
can be represented as M = ±Q z (a)Q x (/3)Q z ( / y), where 
the three independent parameters a, /?, 7 represent the 
three Euler's angles, Qi(x) is a rotation of an angle x 
around the i axis, and the choice of the sign accounts for 
the two possible parity configurations. 

This representation reproduces the symmetric branch- 
ing described in Refs. 12,13 for a = 37r/4, 7 = 7r/4. The 
parameter (3 spans from to tt/2 and represents the cou- 
pling parameter between input and output leads. 

We assume that the effect of the Lorentz force on trans- 
mission and reflection probabilities of the scattering cen- 
ter in Fig. 4b is invariant under cyclic exchange of the 
leads. This symmetry of the system restricts the number 
of free parameters from three to one. Taking 7 as the in- 
dependent variable, the other two angles can be written 
as: 

a = 7r/2 + 7 



sin 27 
2 + sin 27 
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By substitution, we obtain the following scattering ma- 
trix: 

(a b sin 7 b cos 7 
b cos 7 a b sin 7 
6 sin 7 6 cos 7 a 

where the coefficients a and b are given by: 

sin 27 



a± = =F 



2 + sin 27 



6± = ±Vl~a 2 . 

We select solutions a = a + and b = b+ because they 
represent the realistic condition where transmitted elec- 
trons do not gain an extra tt phase. As can be easily seen 
by examining M, the special cases 7 = and 7 = tt/2 de- 
scribe the situation where an incoming electron is totally 
transmitted to one of the outputs or the other. Interme- 
diate cases fall in the range < 7 < 7r/2, the symmetric 
branching condition at zero magnetic field being obtained 
for 7 = 7r/4. 

As a first order approximation, we assume a linear de- 
pendence of the parameter 7 on the magnetic field: 



7 ; 



B 

Bm 



This assumption restricts the allowed values of the 
magnetic field within the range —Bm < B < Bm, where 
Bm represents the magnetic field at which incoming elec- 
trons are fully deflected into only one of the two outputs. 



Electron-phase evolution in the ring arms, including 
electric and magnetic AB phases and decoherence effects, 
is schematized by the black boxes of Fig. 4 and is intro- 
duced in our model by means of 2 x 2 scattering matrices: 
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-* 



here <fi = 2tt^/^o ($0 is the AB quantum of flux h/e 
and the magnetic flux through the ring) represents 
the contribution to the electron-phase evolution due to 
the magnetic AB effect, £ = l/X c the ratio between the 
length of one arm of the ring (I = 2irr/3) and the coher- 
ence length (A c ), and the component due to electron 
propagation along the arm, which amounts to kfl (kf is 
the Fermi wavenumber), and the electric AB effect. In 
our scheme £ models decoherence phenomena occurring 
in our system by decreasing the wavefunction amplitude 
as the electron propagates along the arm. This allows 
to empirically introduce the suppression of the oscilla- 
tion amplitude observed in presence of decoherence at 
the output of QNs in a computationally more efficient 
way compared to the conventional approach consisting 
in introducing random phase shifts with a given proba- 
bility distribution and averaging the output signal over 
it. 

The scattering matrix of the n-th block shown in 
Fig. 4b (n = 1, . . . , 3) composing the ring, including the 
scattering center and the phase evolution due to the prop- 
agation along the arm, is obtained by the combination of 
the two scattering matrices described above: 



J n = I M 2 i e 



•0/3- 



M 12 e 
M 22 e 



M- 



31 



M< 



32 e 



i<j>/3- 
2iG, 
-i(f>/3- 



M< 



23 e 



M13 

i<f>/3- 

NL 



iOn-€ 



33 



We introduced three different phases, # n , to take into 
account that the electron propagation in each arm can 
be independently controlled by means of the three side 
gates: Gi, G2, and G3. In the case where more than 
one one-dimensional subband is available for transport, 
the contribution from each of them can be evaluated by 
calculating n with the corresponding value of kf. This 
would reflect in the total output current as a suppression 
of the AB oscillations, but would not give rise to a qual- 
itative change in the oscillation behavior as a function of 
gate voltage and magnetic field as can be easily verified 
by calculating the transmission and reflection coefficients 
by choosing different values of the parameters n . In the 
following we shall restrict our description to the case of 
a single propagation mode. Since in our experiments we 



keep the gates Gi and G3 at a constant value, we set 
#1 = O3 — 7r/2. As argued above, a different choice would 
have led to a shift in the calculated characteristics of no 
impact in the comparison with the experiment. The re- 
maining phase 62 will be used to model the effect of gate 
G2. The scattering matrix S of the entire ring is calcu- 
lated by combining three of these matrices, one per ring 
node. The probabilities T2 and T3 for an electron to be 
transmitted to the two outputs of the system and R to 
be reflected backwards into the input lead are obtained 
by calculating the squared modules of the components of 
u 7 , given the input condition u = (1,0,0). This situa- 
tion corresponds to injecting electrons in the ring only 
from the input (left) terminal. We wish to emphasize 
that, since in our model electrons that are scattered in- 
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coherently are removed from the system, the sum of the 
reflection and transmission probabilities does not need to 
equal one. 

The impact of decoherence and Lorentz forces on 
the behavior of a three-terminal ring is highlighted in 
Fig. 5, which reports our results for a ring with a ra- 
dius of 160 nm: Panel (a) shows the calculated total 
output probability without these two effects, Panel (b) 
shows the results of the calculations for Bm = 370 mT, 
A c = 320 nm. These values were found by means of a 
best fitting procedure to the experimental data. 
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FIG. 6: Contour plots of the second derivative of the calcu- 
lated transmission probabilities as a function of magnetic-field 
and electric phase shift (O2). (a) Total transmission probabil- 
ity. (b),(c) Individual channel transmission probabilities. 



FIG. 5: Contour plots of the calculated total transmission 
probability as a function of magnetic field and electric phase 
shift (0). (a) Conventional scattering-matrix approach, with 
no Lorentz forces and decoherence (1/Bm = 0, 1/A C = 0), 
(b) Modified scattering-matrix calculations including Lorentz 
force (Bm = 370 mT) and decoherence (A c = 320 nm). 



In Fig. 5a a strong phase-rigidity pattern dominates 
the entire range of external magnetic and electric fields, 
as expected for perfectly coherent, closed systems. On 
the other hand, in Fig. 5b a region reminiscent of phase 
rigidity is present only in a low-magnetic-field region 
(| B I < Bm/3)- At higher field values the oscillation 
phase evolves linearly with the absolute value of B and 
the oscillation amplitude is suppressed due to the Lorentz 
force, as discussed above. It is clear from the comparison 
of the two figures that accounting for the Lorentz force 
and decoherence drastically changes the behavior of the 
total conductance as a function of the external magnetic 
and electric fields. 

In order to allow the comparison of our model results 
to the experimental data, we show in Fig. 6 d 2 T 2 /dB 2 , 
d 2 Ts/dB 2 ^ and their sum, calculated for the same pa- 
rameters of Fig. 5b, as a function of magnetic field and 
gate voltage. The phase of the oscillations reported in 
Figs. 6b and 6c, shows a linear - but opposite in sign 
- dependence from the gate voltage, as observed in our 
devices. 

A further comparison between simulated and experi- 
mental total output of our ring is plotted in Fig. 7. In- 
spection of this figure confirms the ability of our model 



to provide an accurate description of the experimental 
data. 
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FIG. 7: Comparison between the simulated total conductance 
and the second-order derivative with respect to B of the mea- 
sured #2 +#3 as a function of the external magnetic and electric 
fields. 



IV. RANGE OF OPERATION OF A QUANTUM 
NETWORK 



As discussed above, Lorentz force and decoherence 
limit the performance of a QN and should be taken into 
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account when designing such a device. Our theoretical 
model can be exploited to quantify the impact of these 
effects on the operation of a QN. Indeed, by setting the 
minimum acceptable oscillation contrast for the opera- 
tion of a given QN, our model allows one to identify a 
region in parameter space where the chosen requirements 
are met. To this end, it is particularly convenient to fix 
the model variables that are related to external fields 
(parametrized here by (j) and 62) and use as independent, 
dimensionless variables £ and $o/(^mS), where S is the 
ring area, that are related to the sample characteristics. 
Bm can be estimated from the effective arm width W, 
by assuming that almost complete steering of the elec- 
tron in one of the two outputs occurs when r c < W, 
r c = h/(eB) being the cyclotron radius. This condition 
is equivalent to B > and allows to calculate Bm 

as: Bm = &0/W 2 . In the case of the present sample the 
value of Bm yielded by the fitting procedure (370 mT,) 
allows to estimate an effective width W ~ 100 nm, con- 
sistent with the value of the sample arm width (Fig. lb). 
In the following we shall set the minimum visibility at 
15% to match the value of the visibility observed in our 
devices. 

In our model the fraction of electrons that loose phase- 
coherence is removed from the system and can be calcu- 
lated as: 

l-(# + T 2 +T 3 ). 

Also the fraction of coherent electrons that are not useful 
for further computations, owing to the imbalance created 
by Lorentz force and decoherence, must be taken into ac- 
count. This can be evaluated by taking the ratio between 
the amplitude of the AB oscillations and the background 
of the signal transmitted from each output. We evalu- 
ate the contribution of the two effects separately: each of 
them defines a subset of the £, $>o/(BmS) plane, shown 
respectively in Fig. 8a and Fig. 8b. The intersection of 
the two regions (shown in Fig. 8c and calculated at dif- 
ferent values of <£/<£o) corresponds to the set of experi- 
mental parameters that leads to correct operation of the 
QN. Figure 8 shows the operation region taking into ac- 
count decoherence effects (Panel a) and imbalance (Panel 
b). As can be observed from Fig. 8c, requiring the oper- 
ation of the QN within specifications at higher values of 
(fr/^o, results in a more stringent demand on the sample 
parameters, owing to the oscillation visibility decreasing 
at increasing values of the magnetic field. This approach 
can be straightforwardly extended to more complex QNs 
and yields the conditions for a given performance require- 
ment. We believe it can be a useful tool for the design of 
practical QNs. 



V. CONCLUSIONS 

Three-terminal ring devices were fabricated and stud- 
ied at cryogenic temperatures showing clear AB oscil- 
lations as a function of external magnetic and electric 
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FIG. 8: (Color) (a) Subset of the plane identified by £, the 
ratio between the length of one arm of the ring and the co- 
herence length, and by the ratio $o/(BmS), which quantifies 
the impact of the Lorentz force, where the fraction of elec- 
trons that maintain phase coherence when traveling across 
the device is larger than 15%, calculated at a magnetic field 
corresponding to $/$o = 2, and a gate voltage corresponding 
to 6 = 0. (b) Subset of the £, $o/(B M S) plane where the 
fraction of coherent electrons that contributes to the network 
operation is larger than 15%. (c) Subset of the £, &o/(BmS) 
plane where the two conditions depicted in Panels a and b are 
satisfied at = and for different values of $/$o- The filled 
dot represents the values of the parameters £ and &o/(BmS) 
that yielded the best fit to the experimental data. 



fields. The oscillations of the total output signal were 
found to exhibit an evolution from a phase-rigid pattern 
to a linear dependence of the phase on the external elec- 
tric field at increasing (in modulus) magnetic fields. The 
individual outputs displayed a linear dependence of the 
phase on the electric field. We have shown by means 
of a novel scattering-matrix approach that this peculiar 
behavior can be explained only by taking into account 
the competition between decoherence effects and classi- 
cal (Lorentz) forces. Our approach can be successfully 
exploited to determine the range of experimental param- 
eters that guarantee a given oscillation visibility, given 
the geometry and coherence length of a QN. We empha- 
size that this model can be easily applied to more com- 
plex QN systems while still remaining computationally 
efficient. 
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